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We investigate inflation within f{R, (^)-theories, where a dynamical scalar held is coupled to 
gravity. A class of models which can support early-time acceleration with the emerging of an 
effective cosmological constant at high curvature is studied. The dynamics of the held allow for 
exit from inflation leading to the correct amount of inflation in agreement with cosmological data. 
Furthermore, the spectral index and tensor-to-scalar ratio of the models are carefully analyzed. A 
generalization of the theory to incorporate dark matter in the context of mimetic gravity, and further 
extensions of the latter, are also discussed. 

PACS numbers: 


I. INTRODUCTION 

Over the past years, interest towards inflationary cosmology has grown considerably, as a consequence of the great 
amount of data from recent cosmological surveys [in. The inflationary paradigm was first introduced in 1981 by 
Guth 0 and Sato to explain the thermalization of the observable Universe inferred from observations of the 
CMB. It also allows to address some of the problems associated to the initial conditions of a Friedmann universe. 
Moreover, quantum fluctuations during the inflationary epoch presumably seeded the perturbations whi ch g rew 
under gravitational instability into the structures we see today For reviews on inflation, see e.g. [^. 1^. [isl - fTs] . 

An early-time period of acceleration should presumably be supported by a repulsive gravitational force. At the 
same time, a mechanism which allows to quickly exit this stage and enter the radiation dominated era is necessary. 
The arena of inflationary models is quite vast. In the scalar field formulation or chaotic inflation [T^, a scalar field 
(the inflaton) is subject to a potential and drives accelerated expansion when its magnitude is negative and very 
large: at the end of inflation, it settles down in a minimum of the potential and begins oscillating, giving rise to the 
reheating mechanism responsible for particle production. Other implementations of inflation include for instance 
natural inflation (see e.g. Q), k-inflation 0: brane inflation [TTj |. and many others^. In the context of modified 
gravity (see Refs. 0 for reviews), a modification to Einstein’s gravitational action emerges at high curvature and 
supports the early-time acceleration (see Ref. 0). This can be realized for instance in the so-called Starobinsky 
model [l9j |. which provides a correction quadratic in the Ricci scalar. 

A model for inflation is viable only if it is able to reproduce the inferred spectral index and the tensor-to-scalar 
ratio at the origin of cosmological perturbations in the Friedmann universe. The evaluation of these indices depends 
on the theory under investigation. In an unified description has been derived in the context of f{R,cj))-gra.vity, 
where a scalar field subject to a potential is coupled to gravity (a nonminimal coupling in the kinetic part of the field 
is also present). In this work, we will analyze /(i?, (())-inflation by working through some simple examples based on 
modified gravity models which mimic the “false vacuum” of the primordial universe: in fact, we will study a class of 
models (exponential models and power-law models) describing an effective cosmological constant at high curvature 
and whose exit from inflation is induced by the coupling of a scalar field to gravity. We will show that the model can 
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yield values for the spectral index and the tensor-to-scalar ratio (the magnitude of the last one in f{R, (())-gravity is 
particularly small) in agreement with those inferred from observations. 

In the last part of the work we embed this model within the framework of mimetic gravity, which additionally 
endows it with a dark matter candidate. We then discuss extensions which can address the controversies of cold 
dark matter on small scales within the mimetic gravity scenario, such extensions being theoretically driven by 
the correspondence between mimetic gravity and the scalar formulation of the Einstein-aether theories. We then 
speculate on possible further extensions of the scenario depicted. 

The paper is organized as follows. In Section II we will revisit the form of the spectral indices and the tensor-to- 
scalar ratio in /(i?, ^)-gravity by deriving some useful relations. In Section III we study inflation in two different 
/(i?, (/))-models. Early-time acceleration takes place at high curvature in agreement with the latest Planck data and 
the field allows for a quick exit from this stage recovering Friedmann evolution of Einstein’s gravity. In Section IV we 
formulate the model within the mimetic gravity framework. Section V is devoted to conclusions and final remarks. 

We use units where fee = c = /i = 1 and denote the gravitational constant, Gn, by = SttGn, such that 
Gn = Mpi = 1.2 X 10^® GeV being the Planck mass. 


II. INFLATION IN /(i?, <?!))-GRAVITATIONAL MODELS 


Let us consider the following Lagrangian for a scalar field coupled to gravity, 

>- = —-^- 


(II.l) 


where /(i?, </>) is a generic function depending on the Ricci scalar R and the scalar field </> is subject to the potential 
In the above, is in principle a function of (j) which represents a non-minimal coupling of the kinetic term 
of the field. In a flat Friedmann-Robertson-Walker universe, with metric given by 

ds^ = —dt^ -I- a(t)^dx^ , (II-2) 


where a = a{t) is the scale factor depending on the cosmological time t, the Equations of Motion of the theory read: 


3F(i?, + vi<P) + i {RF{R, - f{R, 4,)) - 3HF{R, 4 >), 


(II.3) 


- 2F(R, (t))H = + F{R, 4)) - HF{R, </>). (II.4) 

Here, H = a/a is the Hubble parameter, the dot denoting the time derivative. In the above, we have made use of the 
following definition. 


F{R,4^) = ^f{R,4>). 

From (IH.3I) - (IH.4I) we derive the continuity equation of the field, which reads 


(H.5) 


4> -|- 3F[^ -|- 


1 

2uj{4>) 



dfiR,4>) dV{4>)\ 
d4) d4> } 


0 . 


(H.6) 


Inflation is described by a (quasi) de Sitter solution where the Hubble parameter is nearly constant. To proceed, one 
introduces the “slow-roll” parameters [2l|, 


H 4) ^ F{R,4)) _ E 

iF' “ 2HF{R, 4)) ’ “ 2HE ’ 


(H.7) 


where 


e = e{r,4>M4^) 


?>F{R,4>f 

202 


(H.8) 
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During inflation, the magnitude of the slow roll parameters is presumed to be very small (in what is known as the 
slow-roll approximation). In particular, given that 


- = H + 
a 


(11.9) 


one needs 0 < ei <C 1 in order to have a strong accelerated expansion with H < 0. Correspondingly, the inflation 
epoch and the acceleration end at a time for which ei ~ 1. We observe that 



4:Huj{(j))e3) + 6ei -I- 6£3(1 — £ 2 ) 


0 




HF{R,cl>) ^ 3 



(II.IO) 


where 




—(3 -|- 62 ) H —- 


df{RA) ^dV{4>) 


d(j) 


d4> 


(II.ll) 


Thus, under the slow roll approximation we have that /{H"^F{R, 4>))\ 1 and, since |w((/))^^/F(i?, (/))iJ^| = 

\2uj{(j})^^e3/{HF{R,(j)))\ -C |w((/))(^^/(iJ^F(i?, (/)))|, combining Eq. (|II.3ll with Eq. (|II.6ll yields^ 

3E(i?,0)i/2~E(0) + i(i?E(i?,0)-/(i?,<^)) , (w{cl,)<P - + 2^^^ - 0 . (11.12) 

As a measure of perturbations during inflation, one introduces the spectral index Us and the tensor-to-scalar ratio r 
defined as [^ . 


ris = 1 - 4ei - 2 e 2 -I- 2 e 3 - 2 e 4 , r = 16(ei -I- £ 3 ), 
where £ 1 , 2 , 3,4 must be evaluated during inflation, in the slow roll regime. 


(11.13) 


In the simple scalar field theory with F'{R, (j)) = \/k} and uj{(p) = 1, one has 


ei = 


1 f dV{^) \ 

2k2E(</,)2 \ dcj) J 


£2 = €1 - 


KW((f>) \ d(t> 


d^V{(l)) 


£2 = £4 = 0 . 


Hence, we recover 


rig = 1 — 6£ -I- 277 , r = 16e 


where e , rj are given by 

^ 1 ( dV{<P) y 1 { d^V{cl^) \ 

^ ~ 2KW{(j)Y \ d(l) ) ' \ dcj}^ J ' 


(11.14) 


(11.15) 


(11.16) 


In the modified gravity case instead, with F{R, (j)) = F{R) and w(^) =V{(j)) = 0, given that £2 = 0, £1 ~ —£ 3(1 — £4) 
and, in the slow roll approximation, £1 ~ —£3 and £4 ~ —3£i -I- £i/(i7£i), we obtain 


Us 


1 — 6£i — 2£4 = 1 


2£i 
Hei ’ 


r = 48£ 


2 

1 1 


(11.17) 


where for the expression of the tensor-to-scalar ratio one has to use the higher order corrections of the slow-roll 
parameters. Finally, if = 0 (non-coupling of the field with kinetic energy) and a;((^) = w, the following relation 
holds true 


£1 = -£3 


3HF{R, (j)) 


(£4 -I- 2£3) -I- £3(£2 + £4) — —£3 + 


3HF{R, (j)) 


(£4 -|- 2£3) , 


(11.18) 


^ If Cj{(p) = 0, we obtain directly \uj{(p)^‘^/F{R,(p)H'^\ <C 1 from |e 4 | 'C 1. 
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where we have taken into account the slow roll approximation^. As a consequence, the spectral index and tensor-to- 
scalar ratio read 


Us ~ 1 — 2 ei 


/ 3HF{R, 4>) + 2 a;(^ 2 ' 


- 262 - 6e3 








tlL : 
//2 


0 3F{R, (j)) (HF{R, (j)) - 


^ 2g / 3HF{R,cj))+2oj4'^ \ . 

a;02 j HF{R,<j)) \ 




r = 16(ei + 63 ) = - 


16H 8 F(i?,</)) 


i^2 HF{R,(j))' 


The latest cosmological data from Planck satellite Q constrain these two quantities as rig 
and r < 0.11 (95%CL). 



(11.19) 

( 11 . 20 ) 


= 0.968 ± 0.006 ( 68 % CL) 


III. VIABLE INFLATION IN /(A, (/))-GRAVITY DESCRIBING AN EFFECTIVE COSMOLOGICAL 

CONSTANT 

In order to reproduce the “false vacuum” of inflation, one possibility is to introduce a large effective cosmological 
constant (whose value is close to the Planck scale) within Einstein’s framework. In this way, it is easy to obtain the 
repulsive gravity required to support the early-time acceleration. However, one of the main problems faced by the 
inflationary paradigm is the realization of a mechanism to gracefully exit this stage. 

Over the past years, a class of viable exponential models of /(i?)-modified gravity which can succesfully realize the 
current acceleration of our universe have been investigated. These models feature what can be viewed as a “switching 
on” cosmological constant and assume the following form [23 - 1^ 

/(«) = ^-^( 1 -*■’*') . (m-i) 

where A is the cosmological constant and Rq the curvature scale at which such a constant is expected to appear. It is 
easy to see that for Rq R the model behaves as one where f{R) ~ i? — 2A. Furthermore, by setting A/k^ to being 
the current amount of dark energy in our universe, one recovers the ACDM model. In particular, the dark energy 
epoch is realized by a stable de Sitter solution, with the model passing all cosmological tests. 

It is our intention to extend the model examined in IIII.ll to include inflation. In order to exit from the early-time 
acceleration period (namely, in order to “swith off” the cosmological constant), we will introduce a dynamical field (j) 
by making the following substitution, 


1 

- —bn'^d). 

Ro 


( 111 . 2 ) 


Here, (p is assumed to be negative and dependent on the cosmological time, has been introduced for dimensional 
reasons and 6 is a dimensionless number of order unity. In this way, the scale at which the cosmological constant 
appears is a sort of “running scale”, which varies as the field does. We consider the Lagrangian 


£ = 


i? - 2A (^1 - 



which corresponds to (III.II) with 


f{RA) 


R-2A(^1- 


uj{P) = l, V{P)=0. 


(HL3) 


(111.4) 


® If ^ 0, one has ei ~ —£3 — F{R, <f>)). 
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At the onset of inflation, the field (p = (j)Q is negative and very large, 

(IIL5) 

Given that the curvature during inflation is close to the Planck scale, one may argue that in general 


Mpi < -00 • ( 111 . 6 ) 

The fact that the field is super-Planckian does not pose a problem, provided it is understood that its kinetic energy 
be 0^/2 < Mpi in order to avoid quantum effects in the theory. 

Under the condition set by (IIII.5I) and by assuming the slow roll approximation |ei, 2 , 3 , 4 ! "C 1 in (III.71) . (111.121) allows 
one to get. 


H ~ 



0 ~ 


(111.7) 


where we have used the fact that R = 12H^ + ~ 12H^ and R = 4A becomes the curvature scale at which inflation 

appears. We observe that the field changes faster than the Hubble parameter, as can be noted from (111.41) . where 
H/^ 0^. Thus, inflation ends when 


0 < 


1 

4A6k3 ’ 




R 


( 111 . 8 ) 


and one recovers the Friedmann universe with Einstein’s gravity. In what follows, we will take into account the 
following relation: 

|ij0/(0i7)| ~ l2HH(j)K^ exp [l2H%(t)K^] < 1. (111.9) 

Thus, under the slow roll approximation and by using condition (1111.51) and (1111.91) we get 


F{R, 0) ~ 


F(A,0) ~ 2AAH%'^k^ 


00 + 12b(j)K^ cj)^ 


12H'^b4,K? 


( 111 . 10 ) 


where 

0 ~ . (111.11) 
By combining the expressions in (1111.101) with (1111.71) and (1111.Ill) we obtain 

F{R, 0) ~ 6i7&K^002 , F{R, 0) ~ , 0 ~ 12H%k^^'^ . (111.12) 

Thus, we find that the slow roll paramters in (111.71) are given by: 

ei ~ (4A6 k30o) , £2 ^ IQA %‘^, £3 ~ 16A®&®K^0oe®^^'^°'=", (111.13) 

where we have used (111.41) to evaluate H, and have set i? = 4A and cj) = cpQ. We also obtain 


HF{R, 0) 4A6k®0o 


(111.14) 


As a result, the spectral index and the tensor-to-scalar ratio in (I11.19|) - (lll.20l) with (1111.91) read 

~ 1 - 2£2 = 1 - , r ~ e>(£?_ 3 ). (111.15) 

An important parameter which measures the amount of inflation is the number of e-folds Af, which is defined to be: 




(111.16) 


where ai(U) and af(tf) are the scale factor at the onset and at the end of inflation, respectively, and tjy are the 
respective times. According to the latest data, the number of e-folds must be 55 < A/” < 65 in order for the observable 
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universe to thermalize. Accounting for the fact that the Ricci scalar changes considerably slower than the field itself, 
we can write in our case that 


M = 




Q—4AhK,^(p 

4A6k 


dfj) 




Here, (pi is the value of the field at the end of inflation such that \(pi\ <C |(^o|- Thus, 


l-ris 


2 1 


(111.17) 


(111.18) 


in accordance with the latest Planck data. Since the tensor-to-scalar ratio r is required to be very small but not 
vanishing, by taking into account that the most likely value for r is r ^ 0.06, we require that A^ ~ 10' 


The same analysis may be a ppl ied to another well known f{R) model which mimics the behaviour of (1111.11) . 
namely the Hu-Sawicki model j25| (see [2^ for recent studies of cosmological dynamics within this model), 


fiR) 


R 



1 

l + iR/RoTl ’ 


0 < n. 


(111.19) 


Here, n is a positive fixed parameter and the cosmological constant emerges when 1 <C (i?/i?o)"- To exit from the 
accelerated phase, one can make the substitution (1111.21) in the above expression to obtain 


with the Lagrangian given by: 


f{R,cP) 


R 2A 


1 

1 + {—bcpn^RY 



A 



1 

1 + {—bcpK^R)'^ 



( 111 . 20 ) 


( 111 . 21 ) 


Inflation is realized under the condition given by (1111.51) . namely that the field be negative and its magnitude be larger 
than the Planck scale. Under the slow roll approximation, Eqs. (111.121) can be solved by 


We find 


H ~ 



^ AAHbKn 

{-i2mh(pK^Y+^ ■ 


F{RA) - 
F{RA) - 


2Kb(pKn _ (p^ 

(-12772^,0^3)"+^ ^ 

2M)Kn ( ■■ 12(n + 1)772602,^.3 \ ^ V2{nl)Hb(p^ a ? 

(-1277260a3)"+i + (-127726<).«3) j = {-12 H%(Pk^) 


where we considered the fact that the field varies faster than the Ricci scalar and have made us of 

" 48n(n + l)A77362/t^(^ _ 12(n + 1)77260^^3 

(-127726(?iK3)"+^ “ (-127726 <()k3) 

Thus, the slow roll parameters in (111.71) read 

12n2A262K"‘ 16n(n + 1)A263 k"‘ AA^ii?b^K^ 

'' (-4A60o«3)2(n+i) ’ (-4A6 <)>o«3)"+2 ’ "''^■(-4A6M^)2("+1) ’ 

where (pQ is, as usual, the value of the field at the onset of inflation. We also have 


HF{R,(P) 1 
^2 ~~2' 


( 111 . 22 ) 


(111.23) 


(111.24) 


(111.25) 


(111.26) 
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The spectral index and the tensor-to-scalar ratio in (III.19I) - (III.20I) are found to be 


^ 32n(n + 32 

n,:..!- £2^ ^ = (-4A6</>oAc3)2fo+i) ' 

Moreover, the number of e-folds can be written as: 


M 



4A6Kn 


(—4A6(/)ok^)"'''^ 
16(n + 2)h?h‘^Kf^n ’ 


(III.27) 


(III.28) 


where 4>i is the value of the field at the end of inflation, and we have made use of the fact that the Ricci scalar is 
nearly constant during inflation and changes slower than the field. Thus, 


(1 - Us) = 


2 (n + l) 


8cj)QK 


,2^2 


{n + 2)J\f ’ {n + 2 )^ ' 

In order to satisfy the data from Planck, n must be large. In the limit of 1 <C n one has 


(1 - Us) ~ 


2 ,..2 


A7 


8(/)qK 

JPi 


1 < n, 


and since is on the order of the unit, also the tensor-to-scalar ratio is enough small. 


(III.29) 


(III.30) 


IV. MIMETIC GRAVITY 


Motivated by the recent interest in mimetic gravity, in the present section we aim to contextualize f{R, (/))-gravity 
within such a framework. In Refs. [^ . an approach for expressing the metric in terms of new degrees of freedom 
which isolates the conformal degree of freedom in a covariant way has been proposed (see for further discussions, 
particularly on the role of disformal transformations in mimetic gravity, and related frameworks). The physical metric 
which describes the gravitational systems reads 

gflu = -g^iug°‘^^aip^p^p , (IV.l) 

where is an auxiliary metric and is a scalar field introduced by means of its first derivative. In this way, the 
metric and therefore the action of the theory are invariant under conformal transformations of which take the 
form = Q{t, f2(t, x) being a generic function of the coordinates. As a direct consequence of (IIV.II) one has 

g^‘'df,ipd,^ = -l. (IV.2) 


Such a property can be imposed on the theory by adding a Lagrange multiplier term to the action. Let us consider 
mimetic /(i?, (/))-gravity in (IIIL1|) . with the general action being given in the following form (see also [10, HIl for 
further work on the subject) 


I = 


IM L 






-!/(</.) 


\/-g {gp.u.g^)dsi 


(IV.3) 


In the above, M. is the space-time manifold, while the metric g^^, = g^u^g^v,'^) and its determinant g = g{g^vig^) 
are functions of the auxiliary metric g^i, and the field ip. It should be noted that the auxiliary metric never appears 
explicitly. Here, it is understood that the Ricci scalar is defined with respect to the physical metric g^i, and therefore 
is also a function of of the auxiliary metric and the scalar field if. 


Varying with respect to 5 ^ 1 / and ip yields 

G.. = + T,.) , (IV.4) 

where is the Einstein’s tensor, given by G^i, = — gi_i,jR/2 with the Riemann tensor, and F{R, (j)) is given 

by (|IL5I) . is the stress energy tensor of the scalar field (j) and modifications to Einstein’s gravity are encoded in 
the tensor 


Tt = 


dfj,(t)d^(j) - -gfj,^d°‘(j)da(l) 




(IV.5) 



























(IV.6) 


4>) - F{R, m) + 4>) - g^..aFiR, <^)) . 

Here, and □ are the covariant derivative and the d’Alambertian associated with the metric g^i, = g^iv^ ■ 

Furthermore, the tensor reads 


= - {F{R, </))G -T^- , (IV.7) 

where G , and are the traces of the Einstein’s tensor, the stress energy tensor of (p and of respectively, 

G = -R, 

T* = -Lo{(P)d^(Pd^cP-W{P), 

[2 {f{R, P) - F{R, 4,)R) - 3nF(i?, </.)] . (IV.8) 

By taking the covariant derivative of (IIV.4I) . and recalling that = 0 and = 0, we further obtain that 

[{F{R, <P)G -T*- [y/^ {F{R, p)G -T* - = 0. (IV.9) 

v-g 

Note that the trace of the field equations (IIV.4I) leads to 

{F{R, 4>)G -T^- r“G)(i + = 0, (IV.IO) 

which is automatically satisfied when (IIV.2I) holds true, even when {F{R, (j))G — T'^ — T^'^) ^ 0. 

Here, we stress that it must be 0 < F{R,(j)) to have a positive defined effective Newton constant in = SttGat, 
such that G|f = Gn/F{R,(I)). In this respect, the /(i?, ^)-models (IHI.4I) and (IHI.20I) avoid the “antigravity” in the 
corresponding /(i?)-models (IHI.ll) and (IHL19I) at small curvatures after the end of inflation^ when (/)—>■ 0“ and 
F{R, (p) ~ 1/k^. 

In the context of mimetic gravity, being the field ip not fixed a priori, one is faced with a wider class of solutions, 
as opposed to the simpler case where F{R,p)Gf^u = {TfPi, + A particularly interesting case arises when one 

considers a Friedmann-Robertson-Walker metric (III.21) . which combined with (IIV.2I) yields: 

ip = t + to, (IV.ll) 

In the above, to is an integration constant which can be safely set equal to 0. Hence, Eq.(IV.9) reads: 

- dt {a^ {F{R, p)G - T'^ - t“®)) = 0 . (IV.12) 


From the above we obtain: 


(F{R,P)G-T^-T^^) =-^ , (IV.13) 

where cq is an integration constant. Tf^i, can then be viewed as the stress-energy tensor for an additional component: 

T'm!/= (P + . (IV.14) 

The energy density, pressure and four-velocity of this new component are given respectively by: 

p=-{F{R,p)G-T^-T^^) , p = 0 , u^=d^ip. (IV.15) 

One can then see that T^i, effectively describes a pressureless component, whose energy density decreases as 1/a^. 
The additional degree of freedom provided by the scalar field p, which encodes the conformal mode of gravity, 
can thus mimic cold dark matter. It is dynamical even in the absence of matter (i.e. Cm = 0, as in the case we 


^ Such a problem is not present in the original formulation of these models for the dark energy issue, where the history of the universe 
belongs to Rq < R. 























9 


considered). The amount of dark matter is set by the integration constant cq. 


The mimetic gravity framework offers an alternative approach to solving some of the outstanding problems in 
modern cosmology. On cosmological scales, mimetic dark matter behaves precisely as collisionless cold dark matter, 
and as such is affected by gravitational instability [l^. On the other hand, it is known that the collisionless cold 
dark matter paradigm suffers from a number of shortcomings on small scales: the core-cusp problem [sH, the missing 
satellite problem [s^, and the “too-big-to-fail” problem [s^, just to mention a few (see e.g. [35| for a recent review 
on the subject and more references). While in the particle dark matter framework these issues might be addressed by 
positing that dark matter is self-interacting and collisional (see e.g. [s^ and references therein for further discussions), 
in the mimetic dark matter framework a possible solution is instead the addition of higher derivative (HD) terms for 
the scalar field (p to the action. The HD terms (which are encodin g UV p hysics) provide mimetic dark matter with a 
non-vanishing sound speed, and are effectively dissipation terms^ [37H^ . In fact, these HD terms might be crucial 
to avoid caustic singularities, from which the original mimetic dark matter framework suffers More importantly, 
they have the effect of suppressing power on small scales, which has the potential to address the shortcomings of the 
collisionless cold dark matter framework on galactic and subgalactic scales [13, ■ 

Another important observation is that mimetic gravity is equivalent to a class of Lorentz-violating generally 
covariant extensions of Einstein’s General Relativity, known as Einstein-aether theories (EA hereafter, see [1^ for 
the original formulation of the theory and for more recent extensive discussions on this class of theories). In 
EA, Lorentz invariance is broken by a dynamical unit timelike vector (the “aether”)^ which fixes a preferred 
rest frame at each s pac etime point. In particular, mimetic gravity corresponds to the scalar formulation of the EA 
theory [13, Isee [47| for further discussions), where the aether vector is identified with the gradient of a scalar 
function, = d^ip. This scalar function corresponds to the scalar field encoding the conformal mode in mimetic 
gravity. The most general action for the scalar field, through the inclusion of the aforementioned HD terms, is 
constructed in [45l| . Recently it was also noted that mimetic gravity can be identified and incorporated into the 
framework of covariant renormalizable gravity (48l | . 

Another small-scale open question in mimetic gravity is whether this framework is able to account for the inferred 
flat rotation curves of spiral galaxies (d^ j. The symmetries of the theory allow for the addition of a non-minimal 
coupling between matter and gravity, in the form of a coupling between the aether vector and a matter hydrodynamic 
flux. In the Newtonian limit such a term yields the phenomenology of MOND (see e.g. for comprehensive reviews), 
and hence reproduces flat rotation curves and the Tully-Fisher relation [5l|, among others. The phenomenology and 
constraints on such a coupling remain to be explored. 


The mimetic gravity scenario can be successfully integrated with f{R) gravity, and also with f{R,(j)). Provided 
that the f{R, 4>) (non-mimetic) theory is ghost-free, the corresponding mimetic formulation should also presumably 
be ghost-free, since the addition of a Lagrange multiplier term to the action is not expected to spoil such property 
(although it should be remarked that such a statement remains to be checked) [^. One can then see how 
/(i?, (())-gravity can successfully be incorporated in the mimetic gravity framework to additionally introduce a 
dark matter component to the theory, a crucial element in our current understanding of cosmology. Furthermore, 
it is possible to extend this picture by including a potential for the mimetic scalar field p. Given that in an 
Friedmann-Robertson-Walker universe p can be identified with time, such an addition effectively introduces a 
time-dependent energy density, allowing one to realize any given evolutionary history of the Universe. Any potential 
V{p) can be used to reconstruct a function f{R,(j)) which gives the corresponding evolution, as shown in [30l| . 
In particular, by suitably choosing the form of the potential it is possible to construct an unified and consistent 
description of inflation with graceful exit, the current epoch of acceleration presumably sourced by dark energy, and 
a bouncing non-singular universe. 


In the context of mimetic gravity, our inflationary models (IHI.dl) and (IHI.20I) acquire a dark matter term when 

0 , 


H ~ 



(IV.16) 


® Dark matter with significant dissipation has been studied recently in contexts other than mimetic gravity, both in particle |40| . |41| and 
non-particle frameworks IS respectively. 
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where cq is fixed by (|IV.13I) . The problem with this formulation is that, since the dark matter at the beginning of 
inflation is bounded above by the Planck mass, following inflation its contribution will be completely shifted away 
from the cosmological scenario. Various possible solutions can be envisioned, and here we will briefly discuss two 
possible ones. 

In a first solution, note that when = 0, /(i?, ^)-mimetic gravity leads to Eqs. (III.3ll - pi.4p and we recover the 
same results illustrated in the preceding section. Inflation ends when (^ —>■ 0“ and the models (IIII.4I) and (IIII.20I) 
feature a “phase transition” from a high curvature regime (/(i?, (p) ~ R — 2A) to a low curvature one (/(i?, (p) — R)- 
One may assume that while inflation is realized at = 0, the Friedmann universe belongs to the sector ^ 0 of 
the theory, and dark matter emerges only at the end of inflation. Another possible way of protecting dark matter from 
decay during inflation is to couple the two scalar fields ip and p through a term of the form ipF{(f)), as discussed in (2^. 

Another open question concerns generating the observed radiation and baryonic content in the universe, including 
the observed baryon-antibaryon asymmetry. This can presumably be generated by gravitational particle production 
following the end of inflation, through direct coupling of other fields to the scalar field ip, or through fluctuation- 
dissipation dynamics inherent to the scalar fields (p and ip (see and for further discussions on the topic and 
the implementation of a model of dissipative leptogenesis). We plan to explore these and other ideas in a forthcoming 
paper. 


V. CONCLUSIONS 

In the present paper, we have studied inflation in the context of f{R, (/))-theories of gravity, where a scalar field 
is coupled to gravity. This class of theories is very interesting, given that one can use the /(i?)-gravity sector to 
reproduce a variety of cosmological scenarios (in our specific case, accelerated cosmology at high curvatures and 
Einsteins gravity at low curvatures), while a dynamical scalar field allows for one to move between one scenario and 
another. We note that the /(i?)-formulations of the models under investigation (namely, exponential gravity and the 
so called Hu-Sawiki model with power-law corrections to Einstein’s gravity) belong to a class of viable models for the 
dark energy epoch which the universe undergoes today, where the appearance of an effective cosmological constant 
easily supports an (eternal) accelerated de Sitter expansion. Within the same models (perhaps in the attempt to 
unify the inflationary scenario with the dark energy epoch), one may reproduce the false vacuum of inflation by 
an effective cosmological constant, but a mechanism to make inflation unstable is necessary. In this respect, the 
introduction of a dynamical field induces a phase transition in the models and inflation ends when the effective 
cosmological constant disappears. 

We have explicitly calculated the spectral indices and tensor-to-scalar ratio in the given models, starting from their 
first principle formulation in these kind of theories. We find that for the theory to give the correct amount of inflation 
(namely, a number of e-folds sufficiently large to allow for thermalization of the entire observable universe) and at 
the same time generate a spectral index in agreement with cosmological data, the magnitude of the tensor-to-scalar 
ratio (which is quadratic in one of the slow-roll parameter) is particularly small (note that the same occurs in pure 
modified gravity but not in scalar field inflation within Einstein’s framework) but non-vanishing: this occurs by 
virtue of the fact that the energy scale of inflation is sub-Planckian, while the magnitude of the scalar field can 
exceed the Planck scale. For recent work on /(i?, <())-inflation see also 

The minimal formulation we have considered does not contain a dark matter candidate. To address this point, we 
have then considered extensions of such a model within the mimetic gravity framework, where dark matter appears 
as an integration constant of the equations of motion. In the minimal mimetic gravity formulation, mimetic dark 
matter behaves precisely as collisionless cold dark matter. In the light of the issues which collisionless cold dark 
matter faces on small scales, we have discussed possible extensions of the mimetic gravity framework which allow to 
deal with these shortcomings, and at the same time explain a number of observations, the origin of which is usually 
attributed to particle dark matter (for instance, the inferred flat rotation curves and the Tully-Fisher relation). The 
extensions we have discussed were theoretically motivated by the equivalence between the original formulation of 
mimetic gravity and the Einstein-aether class of Lorentz-violating theories of gravity. We have further expounded 
how the extension of /{R, (p) inflation within the non-minimal mimetic gravity framework allows one to realize 
basically any evolutionary history of the Universe. Finally, we have commented on possible ways of protecting dark 
matter from decay during inflation, and generating the observed baryonic and radiation content of the Universe. 

To conclude, the model we have explored introduces two additional scalar degrees of freedom to the framework of 
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f{R) gravity. A first one allows one to move between two different cosmological scenarios (accelerated expansion and 
Einstein gravity at high and low curvature respectively), while the second one endows the model with a natural dark 
matter candidate (which can address some of the small-scale tensions with collisionless cold dark matter), and can be 
used to reproduce any desired background cosmological expansion. While introducing extra degrees of freedom might 
seem a high price to pay, we have shown that if used appropriately such degrees of freedom allow a more natural 
implementation of a unified expansion history, while at the same time providing a candidate for the missing dark 
matter in the Universe. 
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